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Abstract. Via a non-perturbative method we study the population dynamics 
and photoelectron spectra of Cs atoms subject to intense chirped laser pulses, 
with gaussian beams. We include above threshold ionization spectral peaks. The 
frequency of the laser is near resonance with the 6s- 7p transition. Dominant couplings 
are included exactly, weaker ones accounted for perturbatively. We calculate the 
relevant transition matrix elements, including spin-orbit coupling. The pulse is 
taken to be a hyperbolic secant in time and the chirping a hyperbolic tangent. This 
choice allows the equations of motions for the probability amplitudes to be solved 
analytically as a series expansion in the variable u = (tanh(7ri/T) + l)/2, where t 
is a measure of the pulse length. We find that the chirping changes the ionization 
dynamics and the photoelectron spectra noticeably, especially for longer pulses of 
the order of 10 4 a.u. The peaks shift and change in height, and interference effects 
between the 7p levels are enhanced or diminished according to the amount of chirping 
and its sign. The integrated ionization probability is not strongly affected. 



1. Introduction 

The ionization of an atom by a pulse of laser radiation is a fundamental process which 
is only well understood in certain regions of intensity, frequency and pulse length. The 
regime which is best understood is when perturbation theory may be applied, that is 
where the intensity is very low, the frequency not too low, and the pulse very long 
and constant or quasi constant in electric field amplitude (see e.g. jjj). In particular, 
the amplitude of the initial state is approximately equal to one for the entire duration 
of the process. 

Since the intensity of the radiation field is nowadays often comparable with the 
internal field of the atom, perturbation theory is inadequate. Perturbation theory 
also breaks down at lower intensities if the field frequency is close to a resonance, 
or if saturation occurs. Thus, the applicability of lowest-order theory is limited to 
non-resonant absorption. Resonance effects involving intermediate states can make 
the picture much more complicated and must be described by a theory which allows 
for strong signals. 

Another degree of complexity appears when the length of the pulse is shortened 
(less than 10~ 13 s) : the peaks shift, widen and develop substructures. For example 
the energies of the atomic states shift by amounts, which, to the first non- vanishing 
order, are proportional to the intensity of the laser field (A.C. Stark effect). Thus 
an intermediate atomic level can shift in such a way that it comes into resonance at 
certain times during the laser pulse. 
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The above processes have been approached theoretically in a few different ways 
and described in different, not always compatible languages, often corresponding to 
non identical experimental situations. 

Some calculations are based on the application of Floquet theory ^ |3), a 
semiclassical approach which takes advantage of the temporal periodicity of the 
Hamiltonian. 

Other theories are based on the Keldysh-Faisal-Reiss (KFR) model [Q, |5j. In 
an approximation originally due to Keldysh and developed by Reiss, the exact hnal 
state of the electron is replaced by solutions of the Schrodinger equation for a free 
electron in a laser field. The Keldysh approximation may be quite successful if the 
atomic potential is of very short range, at least for the position of spectral peaks. 
However KFR treatments have not fared particularly well in quantitative comparisons 
with experiments specifically designed to test their predictions regarding ATI electron 
spectra |J. 

Another type of method is the essential states approach which focuses on 
continuum-continuum interactions. States are called "essential" if they are populated 
during the entire process of ATI. Such a theory allows simple analytical solutions, 
provided some assumptions are made. These include the rotating wave approximation, 
the treatment of transitions into the continuum from the ground state as a direct non- 
resonant multiphoton process, and the assumption that all transition matrix elements 
between states belonging to different continua have one of several simple analytic 
forms. 

In a number of studies, the time-dependent Schrodinger equation has been solved 
by direct numerical integration for model atoms. Studies have been done both on 
one-dimensional atoms (e.g. || |[|), and more realistic models using various 
approximations. One of the limitations of this hard numerical approach is the time 
required for these simulations and attendant lack of flexibility. For example the space- 
temporal shape of the laser beam can only be roughly taken into account. This 
approach docs not offer great insight into the physics, but since the experimental 
results can vary greatly even for very similar experiments, these results can be looked 
at as numerical experiments against which one can further test model theories. 

The above-mentioned models usually represent only some important features of 
real systems and make simplifying assumptions. Many studies are done for constant 
fields, or cannot describe processes of resonant-enhanced multiphoton ionization. 
The same can be said about numerical computations, very sophisticated and time 
consuming even for one-dimensional atoms. They only partially contribute to a full 
understanding of the processes. 

Due to the complexities involved in the dynamics of the interaction of an atom 
with strong short laser pulses, it is useful to try to unveil such dynamics using a model 
which is both realistic and offers a simple physical interpretation. 

We report on a model of the system which is solvable, at least partially, by 
analytical methods, since these generally allow greater insight than a direct numerical 
integration. This also reduces the computer effort per parameter set, thus permitting 
the integration over spatial profiles, where many repetitions are needed. 

The work is performed in the semiclassical approximation which has given very 
good results in the past and is justified given the high intensity of the radiation field. 
Furthermore the electric dipole approximation is indeed almost always valid at optical 
frequencies. 

A starting point is the Rosen-Zener problem 111], 03, O . In this problem a two- 
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level atom is subject to a pulse whose amplitude that varies in time as a hyperbolic 
secant. Such problem can be solved analytically via a transformation of the time 
variable t to the "compressed time" u — (tanh(7ri/r) + l)/2. The solutions are 
hypergeometric functions. 

Our model can be regarded as an extension of the basic Rosen-Zener problem to 
an atom with more that two active levels, including coupling to the continuum, to 
model ATI, and chirping of the laser. 

By using an expansion of the wave function in the amplitudes of atomic states 
that are most relevant to the dynamics of the system, one obtains a simplification 
that allows one to get a clearer picture of what is going on during the interaction. 
These levels are those most strongly coupled to the ground state, and others coupled 
to that group. In other words, we use a truncated spectral representation of the wave 
function, which considers only the atomic states most strongly mixed with the ground 
state and includes the contribution of the remaining levels approximately. 

In the application of the method we focus on the cesium atom, both because 
it is advantageous theoretically, since many of the atomic parameters involved have 
already been calculated in the literature, and it is also an element which is accessbile 
to experimental analysis. 

We are interested in studying ATI, thus the coupling of the bound states to the 
continua corresponding to the second ionization peak will also be included. We also 
include frequency chirping in the laser pulse and account for a Gaussian spatial shape. 

Section || introduces the model for the physical system under study and the method 
of solution of the time dependent Schrodinger equation is discussed. In section |^ the 
atomic and laser parameters are characterized. The results of the simulations are 
presented in section [| and the conclusions in section [jj. 

2. Alkali Atoms in Strong Pulsed Resonant Laser Radiation 

In this section we describe the theoretical model and the solutions of the corresponding 
equations of motion. Atomic units are used throughout this article, that is, one takes 
e = m = h = 1 a.u. 

2.1. The Model 

The system under study consists of an alkali atom, which has a single active electron, 
interacting with pulsed laser radiation. The electron wave function is expanded in the 
basis of unperturbed energy eigenstates. To simplify the notation the atomic levels 
are labeled with a single index, which stands for the standard set of quantum numbers 
completely identifying the electron state. In this notation the electron wave function 
can be written as 



where the sum spans the bound states and includes an integration over the continuum. 

As noted in the introduction, we make the dipole approximation. In addition, 
the work is done in the length gauge because it can be argued that it gives a 
better approximation than the velocity gauge when a truncated basis is used to 
represent the electronic wave function. If the unperturbed potentials are velocity- 
independent, the momentum matrix elements are related to the dipole matrix elements 




(1) 



n 
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by < n|p|Q >= i{E n — E ) < n\z\0 >. Thus, the momentum matrix elements decrease 
more slowly with the energy difference than do the length ones. Then, calculating 
multiphoton transition amplitudes, in the sum over intermediate states, the velocity 
gauge weighs states more distant in energy more than the length gauge does and the 
sum converges more slowly. Thus, a truncated basis using states close to resonance is 
a better approximation in the length gauge. The Hamiltonian in length gauge is given 
by 

H = H Q + E-r (2) 

where Hq is the unperturbed atomic Hamiltonian. 
The incident laser field can be written as 

£{t) = 2E(t) cos(flt + a(t)) (3) 

where E represents the pulse shape and a the frequency chirping, i.e. a time- 
dependence of the optical frequency. Experimentally this effect usually accompanies 
the creation of short laser pulses, and to the best of our knowledge, has not been 
included in previous work on ionization. 

We consider photoionization processes in which the laser frequency is nearly 
resonant with the transition between the ground and some excited levels, so that 
ionization is possible from the excited states through the absorption of one photon. 
Specifically we will study cesium subject to laser light in resonance with the 6s-7p 
transition. The level scheme is shown in Fig. [l]. 

If the frequency of the light is close to resonance for the transition between the 
ground state and some excited levels, the dynamics of the electron will be dominated 
by the corresponding contributions of such levels. In addition, other levels should be 
retained if, although nonresonant, they have a very strong coupling with the ground 
state. We call all these states "relevant" states and solve for their effect exactly. 
In addition we can correct for the contribution of the further atomic levels in a 
perturbative way, calculating the energy shifts they induce on the relevant states. We 
use the rotating wave approximation, since we are interested in optical frequencies, a 
regime where such an approximation is typically excellent. 

We make the Weisskopf-Wigner approximation so that the coupling between the 
bound states and the continuum levels is modeled by decay coefficients and level shifts. 
This can be done also for the processes corresponding to excess photon absorption (i.e. 
ATI). 

We now give an outline of the derivation of the equations of motion for the 
probability amplitude using the afore mentioned approximations. Using Eqs.(|l|) and 
(|J) and assuming the electric field to be polarized in the z direction, the time- 
dependent Schrodinger equation yields: 

ic m = £(t) ^ K mn exp(-iuj nm t)c n (t) (4) 

n 

where u) mn = u> m — oj n and K mn =< m\z\n >. Initially the atom is its ground state. 
Since we assume that the non-relevant states and the continuum are going to be 
scarcely populated, we can treat them perturbatively, so that for one of such states, 
k: 

Ck(t) » -i^Kkn / dt'£(t')exp(-w n kt')c n (t'), (5) 
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The sum is over the relevant states. We then substitute the above expression into 
the equation of motion for the relevant state m and perform the time integration 
by parts, we get a term containing a time derivative which we can neglect in the 
Weisskopf-Wigner approximation. 

Further we neglect terms containing factors such as exp(2iQt). We then 
assume that fl w uj nm (quasi-resonant approximation), keep only the resonant 
term explicitly and the first non-vanishing off-resonance correction. Finally, we set 
c m = a m cxp(-zA m J ^ dt'E 2 (t')). 

With these approximations the equations of motion for the probability amplitudes 
{di, i = 1, N} are 

da N 

-jjf = -i^K U E(t)exp[i{A u t + Mt)+a{t))]ai(t), (6) 
i=i 

^ = -i/si ro £?(t)exp[-t(Ai ro t + 0i,(t)+a(t))]ai(t) (7) 
- h m E{t) 2 /2]a m (t), 

<hj{t) = {h-\j) f dt'E(t') 2 (8) 

J — CO 

where the coupling coefficients between levels i and j, i.e. Kij —< i\z\j >, 

Aij = Q — (u)j — uji) is the detuning from the transition frequency, the I summation is 
over all excited states, Xj is the Stark coefficient for the jth level and jjE(t) 2 is the 
decay rate, i.e. "fj = 2ir\njk\ 2 , with uij — u>k = O. Since we will consider only a small 
range of frequencies near resonance, we assume that the decay coefficients are constant 
over this energy interval. The explicit Stark shift terms include contributions from 
all other levels and also the non-resonant contribution from the same level, which 
provides a correction to the rotating wave approximation. For an excited level j, the 
Stark shift is given by 

\ - M 2 | W M 2 | ^ (j + D (9) 



and by 



A _y> \zm\ 2 W / \zik\ 2 \z lk \ 2 \ 



for the ground state, where the sum over fc includes the principal value of the integral 
over the continuum states as well as the bound states not explicitly considered in 
the model, while the sum over n spans the "relevant" states. The equations for the 
continuum amplitude for a state of frequency u K are 

d N 

-J£ = -iE(t) «ifc ex P H + (f>ik{t) + a(t))] a t (t) (11) 
i=i 

for the first ionization peak, and 

, N 

ir - - iE ^ 2 E n« exp H K f + ^ (*) + ai (*) ( 12 ) 



() 



for the second ionization peak (ATI). Here the detuning is given by A' lq = 20 — (ui q — 
u>i), the Stark coefficient = l/f2 2 and xiq is the two-photon coupling between the 
bound state / and the continuum state q (see Eq.(p7|). 

In the above equations we consider only the lowest order terms in the field, i.e. 
cross terms corresponding to the absorption of one photon from an excited state 
and re-emission to a different one are neglected. A few test runs have shown this 
contribution to be negligible. In the same manner, we neglect the effect of the decay 
rate to the second continuum on the equations of motion for the excited bound states. 
One could readily allow for this effect, but it proves to be negligible for the conditions 
we encounter. We also neglect the process involving absorption of two photons from 
an excited state and re-emission of one photon to a continuum state. 

While the equations can be generalized to any number of excess photons absorbed, 
we will consider only the first two peaks. 

To enable the analysis to be carried out in closed form, we chose a field envelope 
in the shape of a hyperbolic secant and a frequency chirp proportional to a hyperbolic 
tangent : 



(13) 



where the characteristic time, r, of the amplitude modulation is large compared to 
This allows a solution in the form of a power series in the compressed time, 
which we will define below. 



E(t) = V4sech(7rt/r), a(t) = (3 / dt' ta,nh(wt'/r) 



2.2. Solution of the Equations of Motion 

In order to solve the equations of motion for the probability amplitude given in the 
previous section, we first let 



bi = ax 

bj = exp[i(Ai m i + 0y(t)+a(i))]a m , 
where m = 2, 3, . . . , N. The equations of motion for the new b amplitude are 

N 

-i^K U E{t)bl{t), 



~dt 



1 = 2 



db m 
dt 



-i Klm E(t)h(t)- 7mE ® b m (t) 
+ i (Ai m + 4> lm (t) + d(t)j b m , 



where 



4>lrn = (A; - \ m )E(t) 2 : 

a(t) = /3tanh(7rf/r). 



(14) 
(15) 



(16) 



(17) 



(18) 
(19) 



The continuum probability amplitudes are found by formally integrating the 
corresponding equations of motion: 



(ik 



(t) = -i / dt'E(t') *lk exp H (Aifct* + <f>i k (t') + a(t'))} k{t')(2Q) 

"'- 00 1=2 
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for the first ionization peak, and 



a q (t) = -i / dt'E(t') 2 J2xk^p[-i(A' lq t' + <f> lk (t') + a(t'))] 6 J (t / )(21) 

J — oo , o 



for the second ionization peak (ATI). 

Letting u = (tanh(7rt/r) + l)/2 and expanding the bound state amplitudes in a 
power series in u, we find : 



(22) 



n=0 



n=0 



n=0 



The expansion coefficients can be calculated through the following recursion equations: 



(j + l)a j+ i + i—Eo V ku/# = 

7T * J 



(23) 



A J) = [n + ^(4A y / + 2/?)- 



n+ I-iJl(A y -/3) 

27^o 2 ^] ^! + (-iAy/ 4^; + 2 7j -I ^)/#2 2 - i^/tyvW (24) 

Using the expansion in powers of u in the formula for the continuum amplitude 
and integrating term by term, we get 



„ oo N 

= -^2^£$>(<W 

n=0 1=2 

x B Q+z^-(A lfc +/3),n+l-z^-(A lfe +/3) 



2r 



x/i U + l- *^-(Aifc - /?), TT + n + i— 2/3, -i(Ai - A fc )/ — (25) 



and 



2tt 



oo N 



2ir 



Z=2 

Z < [ 1 + + 2/3), n + | - ^(Ai fe - 2/3) 



2t 



x 1F1 n + - - z— (Ait - /?), - + n + i— 4/3, -i(Ai - A fe )/ — 02-6) 



2tt 



2tt 



where £> is the beta function and \F\ the confluent hypergeometric function. 
These results have been implemented both with Mathematica and Fortran programs. 
Mathematica allows an arbitrary amount of precision, but it is about a hundred times 
slower than the Fortran version. For fields up to 10 -5 atomic units and pulse lengths 
up to 10 4 atomic units Fortran's double precision proved to yield the same results as 
the Mathematica code. 
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3. Atomic and Laser Parameters 

The model atom is characterized by several coupling and decay constants and level 
shifts, which need to be evaluated. Coefficients for one-photon transitions are given 
by the matrix elements Kfi —< f\z\i > and the Stark-shift parameters A., by Eq.(^|) 
and (JlO|) - The two-photon coupling between the excited states and the continuum is 
given by xji m Eq.(p7|). For the case of the cesium atom the above parameters are 
already available in the literature (see e.g. [|l3)) except for the two-photon couplings 
above threshold. 

We choose the laser field to be linearly polarized in the z direction and of frequency 
such that we can model this problem accurately with the retention of only four excited 
states, those of the 6p and the 7p doublets. The 7p states are quasi resonant, while 
the 6s-6p oscillator strengths are orders of magnitude larger than any other coupling 
to the ground state. Thus, retaining the 6p doublet takes non resonant ionization 
into account with negligible error. The level scheme is shown in Fig. |l|. We take the 
ground state energy to be equal to zero. The ionization potential of cesium is 3.89 eV. 

In Table [j], we list the energy levels and Stark shift coefficients of the relevant 
states of cesium. 



Table 1 . Energy levels and Stark shift coefficients of the relevant states of cesium 



State Energy (a.u.) Stark Shifts (a.u.) 



6>/2 
6 2 Pi/ 2 

6 2 P3/2 
7 2 Pl/2 
7 2 P3/2 



^1=0.0 

^2=0.050931 

o; 3 =0.053456 

1^4=0.0992 

o; 5 =0.1000 



<\ 3 = 
A 4 = 
A 5 = 



-76.93493 

159.221 

232.48688 

64.3317 

386.861 



The decay parameters were calculated in |13| by the method developed by Burgess 
and Seaton |Q, the matrix elements coupling of the bound states were calculated 
using the quantum defect method and the non resonant couplings were calculated using 
data found in Stone pq . In the present work we also need the coupling coefficients 
between the bound states and the continua associated with the second ionization. 
This corresponds to the absorption of 2 photons from the excited states, when the 
absorption of 1 photon is permitted (that is, we have "above threshold ionization"). 

3.1. Calculation of the coupling to ATI channels 
One needs to evaluate terms like : 

Xfi = km V -. (27) 

J , ; ^0 ^ n + UJ, - UJ k + 17] K ' 

The sum spans all intermediate bound and continuum states. Here is the laser 
frequency and Ui is the energy of the initial bound state. The rj in the denominator 
was included because above threshold (fl > \uJi\) the denominator vanishes when 
Lu n — Ui + £1. The small imaginary part prescribes the correct boundary conditions 
for a field adiabatically turned on in the remote past. 

Several methods have been developed to calculate transition amplitudes involving 
infinite sums over intermediate states. For alkali atoms, quantum defect theory and 
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model potentials have been used as weel as truncated summations over a finite number 
of states. Most calculations have taken the laser frequency to be below threshold. 

We used a method based on the solution of an Inhomogeneous Differential Equation 
(IDE), which is equivalent to performing the sums over the entire spectrum of bound 
and continuum intermediate states. In this formulation we have 

X fi =< */|z|*< > (28) 

where > is the solution of 

[H - u]%(r) = -z^i(r), (29) 

with uj = LUi + Q, which is a positive quantity for ATI. 

This method, which is based on numerical computations in the framework of the 
central field approximation has a wider range of applicability than the afore mentioned 
analytical methods, and also leads to accurate results. 

The central potential we used is calculated in ]l6| using the Hartree-Fock-Slater 
approximation. We add a spin-orbit interaction of the form 

V °o(r)=t.*££ (30) 

where V is the Hartree-Fock potential calculated without spin-orbit interaction. We 
work in the SLJM representation. 

The initial state wave function ^i, represents a state identified by the quantum 
numbers i = {n.;, ji,li,m,i} and can be expressed in terms of the product of a radial 
function and a generalized spherical harmonic : 

%(r,£l) = R l {r)y{ Jl ,l u m l -(l)- (31) 

We look for a solution in terms of a superposition of cigcnstates of the total and 
orbital angular momentum. 

%{w,r,tl) = fU(p,j,l,m;r)yV,l,m-,fi). (32) 

j,l,m 

Projecting Eq.( ^9|) on state \jlm >= y(j, I, to; Ci), we obtain : 

[Ho - bj]Ri{ujjlm\r) = - < jlm\ cos9\j l l i m i > R nz {jihmi]r)r. (33) 

The matrix element of cos# can be calculated by an application of the Wigner- 
Eckart theorem Jl7|| : 



(34) 



< l'fm'\r q \ljm > = {-l) 3/2+j +j - m ~ l > ^[{2f + l)(2j + l)] 1/2 r 

r i lUi; r 1/2 

-to' q to J \ j I 1 



where r° = z,^ 1 = =pl/2(x±%) and /> = max(Z,2'). In our case q = 0. 

For cesium the ground state is 6S1/2. Since the final result will not depend on the 
projection of the spin, we can pick a state with ms=l/2 without loss of generality, 
i.e. n= 6, 1=0, j=l/2, m=l/2. Since the interaction z conserves the spin and m;, m 
will always be equal to +1/2 and we will drop it in subsequent formulas. 

After the absorption of one photon from the ground state, the electron will be in 
a state pi/ 2 or P3/2- Since we will consider laser frequencies close to the transition 
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frequency between the ground state and the 7p's states, our atomic model will consider 
these states explicitly, and as noted, also the 6p's, because their coupling to the ground 
state is very large. Thus, the initial state for the transition to the continuum via the 
absorption of 2 photons can be any of the 6p's or 7p's states. Let P(r) = rR(r), then 



dr 2 



1(1 + 1) 



V(jl;r) 



P t {ujl,r) 



< jlh cos 0\jik- > rP n (ji,li;r). 



(35) 



where R n is an eigenstate of the unperturbed atom with n — 6,7, j — 1/2,3/2 
and k = 1. 

First we solve the unperturbed Schrodinger equation to calculate the initial bound 
state which enters in the inhomogeneous term of Eq. (|35|) . This is done through a 
modification of Herman and Skillman's code JT(| to include spin orbit interaction. 

We solve Eq.(|35|) numerically using Numerov's method (see e.g. integrating 
outward from the origin. The boundary conditions require the solution to be zero 
at the origin and to be an outgoing wave for large r. We start the integration using 
a Taylor expansion around the origin. In general the solution obtained does not 
satisfy the boundary conditions at infinity, so we add a solution of the homogeneous 
equation such that the combination satisfies the boundary conditions at infinity. 
The coefficients of the combination are determined using the asymptotic expansion 
of Coulomb functions. In fact, at a large distance from the origin, the potential 
approaches — 1/r so that the solutions have the form of phase-shifted Coulomb waves. 

The transition matrix to a final state \uj q , lf,jf > is given by 

<w f ,l f ,jf\z\y i (w)>=J2< uJf,l f ,j f \z\l'f >< I'j'^iiu) > (36) 
i>ji 

= J2[ drrP u/ (j f ,l r ,r)Pi(uj'l';r) < jfl f \ooa(B)\jl > . 

I'j' J 

The radial integral is calculated by dividing the integration range in two parts. 
The integration is done numerically up to a certain distance i? m ax- The remaining 
integral to infinity is evaluated using an expansion in terms of inverse powers of r, n 
and, whose asymptotic expansion is given in 1 19 as 



Fi{r) = f(r)cos6i - g(r)sm6i, 
Gi{r) = f(r) sinfy + g(r) cos h 
9i =fcr + ln(2fcr)/fc-/7r/2 + cr;, 
ai = argT(Z + 1 - i/k), 



(37) 
(38) 
(39) 
(40) 



where / and g can be expanded in inverse powers of kr for large kr (k = \/2E). The 
resulting integrals have an analytical form, which maybe written as 



dx exp(iax)x x lb 

exp(iai?) \_ ib A — ib | (A — ib)(X — ib — 1) 



-R 



iaR 



[iaRf 



(41) 
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3.2. Transition Amplitudes for Cesium and Laser Parameters 

By virtue of the angular momentum selection rules, the 6s ground state couples to 
the 6p and 7p states, while the p states couple via one-photon absorption to the s 
and d continuum channels and to the p and f continuum channels via two-photon 
absorption. Fig. |^ presents a diagram of the various ionization channels considered in 
this study and Tabs. 0, || and || give the corresponding coefficients calculated with the 
methods described in the previous section. 

Table 2. Coupling coefficients between the ground state and the excited states. 



Parameter 


Value (a.u.) 


«12 


-1.81 


K»13 


-2.56 


ft!4 


0.11 


ft! 5 


0.23 



Table 3. Decay rates from the excited states to the continuum channels via one- 
photon transitions, for a laser frequency quasi resonant with the 6s-7p transition. 
The corresponding coupling coefficients are given by 7ij/27r. 



Parameter 


Value (a.u.) 






721 


117.5878 


(6P1/2 - 


-* S1/2) 


722 


10.45479 


(6P1/2 - 


+ D3/2) 


731 


13.83078 


(6P3/2 - 


■* S1/2) 


732 


10.41901 


(6P3/2 - 


-D3/2) 


733 


93.11069 


(6P3/2 - 


-D3/2) 


741 


37.4357 


(7Pl/2 " 


"* S1/2) 


742 


3.50958 


(7Pi/ 2 - 


-D3/2) 


751 


3.33084 


(7P3/2 - 


-+ S1/2) 


752 


3.003428 


(7P3/2 - 


-D3/2) 


753 


26.8404 


(7P3/2 - 


-D5/2) 



In this work we are interested in short laser pulses. After the first successful 
generation of short optical pulses via mode locking with Ndrglass in the sixties, steady 
progress in generating shorter and shorter pulses was achieved. For a review, see for 
example [pp[ . 

Nowadays the shortest pulses generated directly from a laser can be achieved with 
Tksapphire lasers. These lasers can be used to achieve 10~ 14 s pulse generation |H| 
by using nonlinear effects which induce frequency chirps. Complete compensation of 
the chirp is not possible with a finite number of linear optical elements. 

With these experimental considerations in mind, we will consider laser pulses 
of durations in the 25-250 fs range, with frequency almost on resonance with the 
6s-7p transition. The frequency chirping induced by the compression mechanisms, 
is characterized by a chirp parameter [3 of the order of tt/t. The parameters 
characterizing the laser beam are thus frequency, chirp parameter and intensity. The 
temporal dependence of the pulse is given by Eqs.(2.3) and (2.16). We will consider 
intensities in the range 10~ 7 — 10~ 4 a.u., and pulse lengths of 10 3 — 10 4 a.u., with chirp 
parameters /? w tt/t. 
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Table 4. Coupling parameters to the continuum via two-photon transitions (ATI 
peak); i is the imaginary unit. 



Parameter Value (a.u.) 



X24 


1972 - 2916i 


(6P1/2 - 


-P1/2) 


X25 


-960 + 484i 


(6P1/2 - 


-P3/2) 


X26 


-2568 + 4759i 


(6P1/2 - 


-F5/2) 


X34 


-616 - 884i 


(6P3/2 - 


-P1/2) 


X35 


3188 + 2068i 


(6P3/2 - 


- P3/2) 


X36 


844 + 456i 


(6P3/2 - 


^F 5/2 ) 


X37 


-2672 - 1356i 


(6P3/2 - 


-F 7/2 ) 


X44 


106.8 + 29i 


(7Pi/ 2 " 


-P1/2) 


X45 


64 - 15.9i 


(7Pi/ 2 - 


-P3/2) 


X46 


383.2 - 460i 


(7Pi/ 2 - 


-F 5/2 ) 


X54 


62.8 - 26.84i 


(7P3/2 - 


-P1/2) 


X55 


-142.8 + 23.92i 


(7P 3 /2 - 


-P3/2) 


X56 


-65.6 + 79.6i 


(7P3/2 - 


- F 5/2 ) 


X57 


199.2 - 235. 6i 


(7P3/2 - 


^F 7/2 ) 



4. Ionization Dynamics of Cesium 

In this section, the method described in section is applied to the case of cesium, 
using the atomic and laser parameters calculated and listed in section ||. We will 
consider three different regimes according to the intensity of the laser beam and the 
case of Gaussian beams. The 'atomic' unit of intensity here is 1.4 x 10 17 W/cm 2 - 
four times the unit used by many other authors. 

4-.1. Weak Field Limit 

The spectrum for weak fields (in the adiabatic limit) is usually a single peak at an 
energy twice the photon energy above the ground state. In fact, for this case, one 
expects the final energy of the photoelectron to consist of the Fourier distribution of 
the pulse squared at an energy consistent with energy conservation. However, with 
relatively short, weak pulses tuned between the doublet, a double peak structure is 
found in the Si/ 2 or d 3 / 2 continuum channels but not in the d 5 / 2 channel. This 
structure was first reported in fi~3| and attributed to the interference between the 
two different ionization paths to the afore mentioned channels, which results in the 
suppression of the central peak. However the d 5 / 2 channels can be reached only in one 
way, through the p 3 / 2 excited states and thus cannot show any interference effects. 

Adler et al. jl3| found that one sees a pair of peaks only when the condition 
1 < |Aijt| < 10 is roughly true. This behavior is consistent with the interpretation 
of such structure as a manifestation of destructive interference between the 7p levels, 
when the laser is tuned to the middle of the doublet. For very short pulses, the 
incident radiation's spectrum covers a wider range of frequencies and the detuning 
is less significant so that the 7p's are practically a new single state and the doublet 
disappears (see Fig. ||). For longer pulses, the central peak forms again and hides the 
side peaks, which are much smaller. The central peak is about at frequency uj\ + 2 ft, 
while the two side peaks at U4 + £1 and + f2. The relative heights are correlated 
with the average populations of the 7P levels during the pulse. 

In Fig. U we show the spectra of the Si/ 2 channel and the the spectra corresponding 
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to the contribution from the 7p states separately, i.e. with no inteference. One can 
see here that each 7p state contributes to one of the side peaks. The central peak is 
not resolved, being very close to the side peaks. If there were no interference effects 
we would not have the dip between the two side peaks. The existence of two peaks is 
not due do oscillations in the populations, since we find that these do not occur for 
the intensities we are considering here, as can be seen for example in Fig. ^. 

In the ATI peak the double peak structure is usually absent, in fact the mediation 
of several channels to the first continuum, makes interference effects less likely to 
occur. For some particular frequencies there is still some structure present in the P3/2 
channel, which is the channel receiving contributions from all first peak channels. See 
Fig. ||, where we plot ATI spectra for the P channels for different pulse lengths. As 
for the two-photon peaks, also here the structures are washed out for shorter pulses. 
In general the ATI peaks are at frequency uo\ + 3£! and have a bandwidth larger than 
that of the lo\ + 2fl peaks. Infact they entail the Fourier transform of the field pulse 
to the third power, which has a smaller width (in time) than the square of the pulse. 

We note that the broadening of the ATI peaks was reported in the experimental 
results reported in 1992 by Nicklich et al in f22| . The calculation described by these 
authors did not predict broadening of the higher peaks, however, although their spectra 
did contain the oscillations similar to those in the experiment, while those oscillations 
are absent from our work. 

The double peak structure is heavily influenced by the presence of chirping in the 
laser light, especially for the longer pulses in the above range. Varying the strength 
and the sign of the chirping changes the absolute and relative height of the peaks, see 
Fig. |j. In fact, the height of the peaks is determined by the amplitudes of the excited 
states and their phase relationship, which are influenced by a kind of self-induced 
phase modulation. Thus imposing external phase modulation (through chirping) on 
these can strongly influence the ionization probability. 

The evolution of the population of the bound states also changes appreciably with 
chirping for the 7p levels which arc quasi resonant with the laser radiation. The 6p 
levels are only slightly affected (see Fig. ^). Infact the 6p's are off-resonance and 
their populations follow the field adiabatically. In addition the relative height of the 
side peaks is related to the amount of population in the corresponding level. The sign 
of the chirping determines which levels will have a bigger average population during 
the pulse and thus the height of the spectral peak. With our conventions a positive 
chirp corresponds to a pulse whose instantaneous frequency is lower than the carrier 
frequency at the beginning of the pulse and increases to a value greater than the carrier 
frequency in the second half of the pulse. Thus the 7pi/2 is closer to resonance at the 
beginning of the pulse and starts becoming populated earlier in the pulse. The 7p3/2 
on the contrary becomes populated at a higher rate later during the pulse, when it 
gets closer to resonance. Thus its average population during the pulse is lower than 
in the case of chirp free light. The situation is reversed for negative chirping. 

4- 2. Strong Field Limit 

In the previous section we observed that the double peak structure disappears for 
short pulses. However, Adler et al. ^3| have shown that when the pulse is short and 
the intensity of the beam increased, the spectra again contain a multipeak structure, 
which is also present when the laser frequency is tuned outside the doublet. Multiple 
peaks can appear also in the D 5 / 2 channel. We can see this in Fig. ||, where we plot 
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the channel D 3 / 2 spectra for different intensities. Increasing the laser intensity we see 
a multipeak structure emerge, at the same time the populations also start showing 
oscillations, see Fig. ||. 

As noted in the previous subsection, this result confirms that other mechanisms 
such as Rabi oscillations and Stark shifts play a role in the production of the multipeak 
structure in strong fields. 

As in the case of weak fields, for such short pulses, the influence of chirping 
is minimal in the second continuum peaks, but it induces visible changes in the 
population dynamics and the relative peak heights of the doublet. 

For long pulses, the chirping causes dramatic changes in the population dynamics 
and in the spectra of the threshold peaks, see Figs. [r| H The ATI peaks are instead 
not strongly affected as one expects because of the larger linewidth and general lack 
of structure, as shown in Fig. |l2|. It appears that chirping enhances or quenches the 
oscillations in the populations depending if it is an up-chirp or a down-chirp, while 
the frequency of the oscillation remains the same. 

4-3. Gaussian Beams 

In actual experiments the laser beam incident upon the atomic sample has an intensity 
spatial profile, so that the peak field seen by each bound system is a function of its 
distance from the axis. In weak fields this is not important, since the intensity is 
merely an overall factor. This is not the case in strong fields, however. Both the 
total ionization and the electron energy distribution become complicated functions of 
the optical field strength, as for example, in Fig. || in which we see a second peak 
developing in the D 3 / 2 channel with increasing intensity. 

If the intensity of the field is weak enough to allow the neglect of depletion of the 
ground state and coupling between the excited states through multiple transition to 
the ground state and back, then none excited states is affected by the existence of the 
other ones. In this limit the solution of each bound state amplitude takes the form 
of the Rosen-Zener solution (see for example |l^] and jl]J ) for a two-level system in 
which all energy shifts have been neglected. The spectra produced by this expression 
are identical to those produced by the exact result for weak laser intensities. 

We examine for what intensities and pulse lengths perturbation theory breaks 
down. The evolution of the bound states start to deviate from the typical behavior 
at weak fields when the peak intensity and pulse length are such that J max T w 10 
a.u., for which nonlinear effects start to show. This onset of nonlinear effects seems 
to be approximately independent of chirping. 

In the perturbative regime, the total ionization is proportional to the square of 
the intensity, since the dominant process is two-photon absorption. 

The total ionization probability is not very sensitive to the presence of chirping, 
(Fig. |r|). There is in general a small increase in the ionization rate, which can be 
attributed to the widening of the pulse bandwidth from the frequency chirp. The 
scaling parameter is 27 max T/7r and the total ionization is practically the same when 
the length of the pulse or the current / is changed maintaining the scaling parameter 
constant. 

In the figures we show the evolutions of the populations and the spectra as a 
function of the distance from the center of the laser beam. 

Figs. [l4| show the total ionization probability and the excited state population 
after the passage of the pulse, as a function of the distance from the center of the 
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beam. A curve showing the Gaussian shape of the beam spatial profile is shown for 
reference. We find that in general the total ionization probability roughly follows the 
intensity profile of the laser beam, while the population of the excited states after 
the passage of the pulse shows an oscillatory dependence on the laser peak intensity 
at the location of the atom. Looking at the temporal evolution of the excited states 
populations one can see that the gradual change is due to the change in the Rabi 
frequency and the number of oscillations the populations can do during the pulse. As 
an illustration in Figs. 15 and we show the population of the 7p 3 / 2 state as a function 
of the distance from the center of the laser beam without chirping and with chirping. 
The distance is normalized such that the FWHM of the beam is unity with the laser 
parameters given in the figure caption. The effect of chirping is mainly to change the 
height of the oscillations, a chirping of the opposite sign would flatten out the surface 
of the plot. In Fig. [TJ] and we show the spectra for the sj/2 channel for the same 
set of parameters as Figs. |l5| respectively. The presence of chirping here changes the 
shape and heights of the peaks. One can see that for this intensity and frequency the 
spectra does not depend in a monotonic way on the intensity. The final spectra will 
be a spatial average of all these contributions and the final population of the bound 
states will also change with position. Then in Fig. ^ we show a comparison between 
the spatially averaged spectra and the spectra at the peak intensity. 



5. Conclusion 



The goal of this work was to model the dynamics of an atom with a finite number of 
active levels interacting with short pulse of strong laser light using a non-perturbative 
approach in the dominant states. It was hoped to strike a balance between the 
competing needs of realism and numerical simplicity, extending prior analysis to the 
case of ATI and a "chirped" pulse. 

As was done in prior work, we assumed only a few active levels whose coupling 
to other states was taken into account by means of time-dependent level shifts and 
decay rates. By assuming a hyperbolic secant pulse shape for the amplitude of the 
field, and a related form for the frequency modulation, we were able to extend the 
analytically solvable problem previously reported to the present case. The ATI was 
accounted for by means of an effective operator, whose coupling strength was obtained 
from the solution of an in- homogeneous Schrodinger equation. 

The model was applied to the cesium atom for light nearly resonant with transitions 
between the 6s ground state and 7p excited doublet. Coupling to the 6p doublet is 
also included explicitly, since, even though its frequency is far from resonance, the 
interaction matrix element is very large. The approach could be further extended 
to include more atomic levels and a larger number of ATI peaks. Allowance for the 
temporal profile of the laser beam was made, and the dynamics of the excited state 
populations and spectra of the ejected electrons were efficiently calculated. 

We found that chirping as a great influence on the dynamics of the quasi- resonant 
states, but much less of an effect on the off resonant states. This transpires because the 
chirping causes much less of a fractional change in the detuning of the nonresonant 
states. In fact, it appears that the dominant part of the change produced in the 
nonresonant amplitudes by the chirping is due to the modification in the ground state 
amplitude cause by the resonant intermediate couplings. 

The effect of chirping is more pronounced for longer pulses. In this case, the pulse 
has a smaller linewidth the the fractional change due to the chirping is relatively more 
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prominent. 

Some of the spectral features, such as double peaks, can change drastically with 
the chirpings. Those effects in weak fields are due to the interference between the 
phases of the various atomic amplitudes when the laser is tuned between resonances, 
and is very sensitive to time dependent frequency shifrts. The interference is enhanced 
or suppressed depending on the sign and magnitude of the chirp. 

The structure of the spectral peaks is quite sensitive to the length of the pulse. 
The doublet structure disappears for short pulses. The shorter pulses have a broader 
frequency spectrum, over which the populations dynamics are averaged, thus washing 
out the deleicate interference effect due to phase modulation. 

It is interesting to note that the ATI peaks lack the structure of the principal two- 
photon peaks, even though both types are generated through the same quasi-resonant 
intermediate states. This appears to result from the following considerations. The 
weak-field structure in the the two photon case stems from the interference between 
the pl/2 and p3/2 amplitudes in a case when they are tending to cancel, so that 
the overall transition probability is small. In the three-photon peak, there are more 
channels, so that the net probability is dominated by the contributions that are not 
supressed by destructive in- terference. 

However, according to our calculations, the three-photon peak is broader than 
the two photon by about a factor two. This can be understood in terms of the 
Fourier transforms of the effective operator coupling the intermediate states to the 
respective continua. For two photon ionization, that operator is just the transform of 
the hyperbolic secant, while for the ATI peak, it is the transform of the square of the 
hyperbolic secant, which is broader in frequency space. 

We note that chirping affects not only the shape of the photoelectron structure , 
but also the total ionization probability, perhaps by making more frquencies available 
for ionization. 

We also studied the effect of spatial dependence of the laser beam profile on the 
quantities accessible by experiment. The total ionization, the spectra and the residual 
population show deviations from perturbative behaviour which should be taken into 
account performing a spatial average, to find the measured quantity. The effects 
of a gaussian beam profile was discussed in J|4| for the case of a one-dimensional 
atomic model. There the concurrence of significant amounts of ionization and 
residual excited-state populations is attributed to the fact that atoms are subjected 
to different peak intensities depending on position, also oscillations in the residual 
excited populations are found as a function of intensity. 

As discussed previously we also find that for a certain range of intensities and 
pulse lengths, the total ionization is proportional to the product of peak intensity 
and pulse length. This seems to imply that the total ionization is behaving like a 
one-photon absorption out of the p states, and that this step is roughly given by first 
order perturbation theory. 

Assuming a Gaussian beam, the space averaged spectra can show qualitatively 
different features from the spectra corresponding to the peak intensity. In fact some 
of the spectral features are very sensitive to the peak intensity. The residual excited 
state population shows an oscillatory dependence on the laser intensity, while the total 
ionization increases monotonically with the intensity. 

In summary, we have shown that the methods used by Adler et al. can be extended 
to the case of chirped pulses and to the inclusion of ATI. The population dynamics 
and ionization spectra of cesium atoms subjected to strong short laser pulses has been 
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analyzed. We have found that the behavior of that system is quite sensitive to the 
laser light characteristics, such as pulse lengths, frequency chirp and intensity spatial 
profile, that are not usually taken into account. 
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Figure 1. Energy level structure for the Cesium five level system, with the excited 
states coupled to a continuum of states 




Figure 2. Ionization channels from the ground state of cesium 
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Figure 3. Probability per unit energy vs. energy for the photoelectron produced by 
a two photon ionization of cesium for various pulse lengths. Laser frequency 0.0994 
a.u., peak intensity 10~ 7 a.u. (1.4 X 10 10 W/cm 2 ), no chirping. 




Figure 4. Probability per unit energy vs. energy for the photoelectron produced 
by a two photon ionization of cesium. The solid line is the spectrum for the S1/2 
channel, the other curves represent the contribution of ionization from the 7p levels 
separately (no interference effects). Laser frequency 0.0994 a.u., peak intensity 
10~ 7 a.u. (1.4 X 10 10 W/cm 2 ), pulse length 250fs, no chirping. 
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Figure 5. Probability per unit energy vs. energy for the photoelectron produced by 
a two photon ionization of cesium for various pulse lengths. Laser frequency 0.0997 
a.u., peak intensity 10~ 7 a.u. (1.4 X 10 10 W/cm 2 ), no chirping; pulse length r is — 
100 fs, 150 fs, ■■ 200 fs,— 250 fs. 




Figure 6. Probability per unit energy vs. energy for the photoelectron produced 
by a two photon ionization of cesium for different chirping parameters. The final 
continuum state is an S1/2 channel. Laser frequency 0.0994 a.u., pulse length 250 fs, 
peak intensity 10~ 7 a.u. (1.4 X 10 10 W/cm 2 ), line type correspond to different chirp 

parameters 0: - -0.0003 a.u., 0.0002 a.u., o 0.0 a.u., •• 0.0002 a.u., -■ 0.0003 

a.u. 
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Figure 7. Populations of the 7p levels vs compressed time for different chirp 
parameters. Laser frequency 0.0994 a.u., pulse length 250 fs, peak intensity f0~ 7 
a.u. (f.4 X f0 10 W/cm 2 ), line type correspond to different chirp parameters (3: — 
-0.0003 a.u., 0.0002 a.u., o 0.0 a.u., •• 0.0002 a.u., -• 0.0003 a.u. 
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Figure 8. Probability per unit energy vs. energy for the photoelectron produced by 
a two photon ionization of cesium for various intensities. The final continuum state 
is an D3/2 channel. The laser is tuned outside the 7p doublet, at the frequency 0.099 
a.u., the pulse length is 25 fs. The peak intensity in the legend is given in atomic 
units. 




Figure 9. Population of the 7p levels vs compressed time, for various intensities. 
The parameters are the same as in fig. M 
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Figure 10. Populations of the 7p levels vs compressed time for different chirp 
parameters. Laser frequency 0.0994 a.u., pulse length 250 fs, peak intensity 

0.7 X 10~ 5 a.u., line type correspond to different chirp parameters /3: 0.0003 a.u., 

0.0002 a.u., o 0.0 a.u., ■• 0.0002 a.u., -• 0.0003 a.u. 
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Figure 11. Probability per unit energy vs. energy for the photoelcctron produced 
by a two photon ionization of cesium for different chirping parameters. The final 
continuum states are the S1/2, D 3 / 2 and D5/2 channels and the total probability for 
2 photon ionization. Laser frequency 0.0994 a.u., pulse length 250 fs, peak intensity 

0.7 X 10~ 5 a.u., line type correspond to different chirp parameters f3: 0.0003 a.u., 

0.0002 a.u., solid/dotted line 0.0 a.u., ■■ 0.0002 a.u., -■ 0.0003 a.u. 
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Figure 12. Probability per unit energy vs. energy for the photoelectron produced by 
a three photon ionization of cesium for different chirping parameters. Laser frequency 
0.0994 a.u., pulse length 250 fs, peak intensity 0.7 X 10 — 5 a.u., line type correspond 

to different chirp papramctcrs (3: 0.0003 a.u., 0.0002 a.u., solid/dotted line 

a.u., ■■ 0.0002 a.u., -■ 0.0003 a.u. 
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Figure 13. Total ionization probability as a function of intensity for different 
chirpings (in a.u.). The laser frequency is 0.0995 a.u. and the pulse length 250fs. 
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Figure 14. Residual total population of the excited states and total ionization 
probability vs. distance r from the axis of a Gaussian laser beam. The Gaussian 
beam profile is shown for reference (dash-dotted line). Laser frequency 0.0997 a.u., 
pulse length 250 fs, peak intensity 0.7xl0 12 W/cm 2 , no chirping. 




Figure 15. Population vs. compressed time u and distance r from the axis of a 
Gaussian laser beam. Laser frequency 0.0993 a.u., pulse length 250 fs, peak intensity 
0.7xl0 12 W/cm 2 , no chirping. 
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Figure 16. Probability per unit energy vs. energy for the photoelectron produced 
by a two photon ionization of cesium vs. energy Ef and distance r from the axis of a 
Gaussian laser beam. Laser frequency 0.0993 a.u., pulse length 250 fs, peak intensity 
0.7xl0 12 W/cm 2 , no chirping. 



S1/2 channel 




196 0.197 0.198 0.199 0.2 0.201 0.202 

Final Energy 

Figure 17. Probability per unit energy for the photoelcctron produced by a two 
photon ionization of cesium vs. energy Ef. The final continuum states are the S1/2 
Laser frequency 0.0993 a.u., peak intensity 0.7 xlO 12 W/cm 2 , pulse length 250 fs, 
no chirping. Space averaged spectra (continuous line) and peak intensity spectra 
(dashed line). 



